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Abstract. We consider the effect of energy dependent 
propagation of relativistic electrons in a spatially inho- 
mogeneous medium in order to interpret the broad-band 
nonthermal radiation of the young shell-type supernova 
remnant (SNR) Cassiopeia A. A two-zone model is pro- 
posed that distinguishes between compact, bright steep- 
spectrum radio knots and the bright fragmented radio ring 
on the one hand, and the remainder of the shell - the dif- 
fuse 'plateau' - on the other hand. In the framework of 
this model it is possible to explain the basic features of 
the spectral and temporal evolution of the synchrotron 
radiation of Cas A if one assumes that these compact 
structures correspond to sites of efficient electron acceler- 
ation producing hard spectra of accelerated particles with 
power-law indices /3 acc ~ 2.2. The resulting energy distri- 
bution of radio electrons in these compact structures be- 
comes significantly steeper than the electron production 
spectrum on timescales of the energy dependent escape 
of these electrons into the surrounding diffuse plateau re- 
gion. We argue that the steepness, rather than the hard- 
ness, of the radio spectra of compact bright structures in 
clumpy sources can in general be considered as a typical 
signature of sites where strong electron acceleration has 
built up high gradients in the spatial distribution of radio 
electrons. Subsequent diffusive escape then modifies their 
energy distribution, leading to potentially observable spa- 
tial variations of spectral indices within the radio source. 
Qualitative and quantitative interpretations of a number 
of observational data of Cas A are given. Predictions fol- 
lowing from the model are discussed. 

Key words: acceleration of particles - radiation mecha- 
nisms: non-thermal - radiative transfer - supernovae: in- 
dividual: Cas A - radio continuum: stars 



1. Introduction 

Cassiopeia A is the youngest of the known galactic super- 
nova remnants (SNRs) whose birth probably dates back 



to 1680 (Ashworth 1980). It is one of the most prominent 
and w ell st udied radio sources on th e sky (e.g. Bell et al. 
1975| ; Tuffs |l986| ; Braun et al. |l987|; An derson et al. |l991 



- hereafter ARLPB; Kassim et al. |1995| , hereafter KPDE; 
etc.), whose synchrotron radiation pr obabl y extends into 
the h ard X-ray region (Allen et al. 1997 ; Favata et al. 
1997 ). Most of the radiation of both nonthermal and ther- 



mal origin comes from a shell region enclosed between two 
spheres, with angular radii ~ 150 arcsec and ~ 100 arcsec, 
corresponding to spatial radii Rq = 2.5 pc and i? r i ng = 
1.7 pc, r espectively, for a distance of d ps 3.4 kpc (Reed et 
al. 1995). The former corresponds to the mean radius of 



the assumed blast wave, while the latter is supposedly the 
mean radius of the reverse shock in the freely expanding 
(upstream) ejecta which are heating the gas to tempera- 
tures ~ (1 — 3) keV, thus creating the hot thermal X-ray 
component (Becker et al. 1979; Fabian et al. 198C; Jansen 
et al. |i~98Sl). 



Cas A is considered to be a powerful particle acceler- 
ator, with an energy content in relativistic electrons, esti- 
mated from simple equipartition arguments, in the range 
Wo - 10 48 - 5 x 10 49 erg (e.g. Chevalier et al. |1978[ ARLPB 
) . This value should be increased by more than an order of 
magnitude for the total energetics in relativistic particles 
if the high ratio between relativistic protons and electrons 
observed in cosmic rays (CR) holds in Cas A. However, 
the mechanisms and sites of particle acceleration in Cas A 
have not yet been identified. 

Observations of Cas A at optical wavelengths show 
strong line emission from 2 main types of compact 
structures. The quasi stationary flocculi are thou ght t o 
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trace the circumstellar medium (e.g. Fesen et al. 1988 ), 
while numerous fast moving knots (FMK) represent dense 
clumps of supernova ejecta moving ballistically with aver- 
age velocities of ~ 53 00 km /s from the time of explosion 
(see e.g. Reed et al. 1995| , and references therein). The 
FMKs have been invoked as sites of particle acceleration 
by Scott & Chevalier ( 1975| ) who proposed second order 
Fermi acceleration of particles in the turbulence created 
in the wake of FMKs overtaking the shell. Another pos- 
sibility is a first order Fermi acceleration process at the 
bow shocks driven ahead of FMKs as proposed by Jones 
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et al. ( 1994 ). These authors showed that these dense gas 
'bullets' can effectively accelerate electrons at the stage of 
deceleration and subsequent fast destruction of the FMKs 
by Kelvin-Helmholtz and Rayleigh- Taylor instabilities. 

Perhaps the most straightforward suggestion, however, 
is that efficient electron acceleration occurs directly in re- 
gions of high radio brightness - in other words, that these 
regions are bright not only because of enhanced magnetic 
field, but also due to local enhancement of relativistic elec- 
trons. Such regions have been suggested to be sites of elec- 
tron acceleration by a number of authors (e.g. Bell 1977, 



Dickel & Greisen |1979J Cowsik & Sarkar [198J ). In Cas A 
they are represented by rather compact structures which 



include > 300 radio knots (Tuffs [L98§ ARLPB) located in 
the shell, several paraboloidal radio features suggested to 
be bo w sho cks associated with decelerated ejecta (Braun 
et al. 1987), and the bright fragmented radio ring at the 



projected radius of the (supposedly) reverse shock/contact 
discontinuity. Examples of all these morphological struc- 
tures are depicted in Fig.l. 

Any theory of particle acceleration in Cas A has to 
address a wide range of spatial and spectral characteristics 
of the observed nonthermal radiation from the radio to 
infrared (IR), and possibly to X-ray regimes, which can 
be summarized as follows: 

— For 20 MHz < v < 30 GHz the total radio flux is well 
approximated by a power-law J{v) oc v~ a with an 
index a w 0.77 (e.g. Baars et al. |1977[ ); 

— For v < 20 MHz the spectrum J(v) essentially flattens, 
and turns over below 15 MHz (Baars et al. 1977); 

— The secular decline of the fluxes seems to be frequency- 
dependent, with a decline rate at the level of (0.8 — 
1) %/yr in the frequency range v ~ 40 MHz — 1 GHz 
(Rees |1990[ Hook et al. 1 992), and - 0.6 %/yr at v ~ 
10 GHz (Dent et al. [l974] O'Sullivan & Green ^999| ). 

— The contribution of the diffuse emission of the shell 
(plateau) to the total flux of Cas A at 5 GHz amounts 
to w 50 %, and the second half is due to the frag- 
ments of bright radio rin g (30 %) and radio knots + 
bow shocks (20%) (Tuffs |l 98 6| ). The total flux at the 
epoch 1987 was about 750 Jy (see ARLPB). 

— Individual radio knots show a wide spread of spectral 
indices in t he ran ge fro m a ~ 0.6 to a ~ (0.9 — 0.95) 
(Rosenberg [l97C| ; Tuffs |l983j ; ARLPB; AR96). 

— For v > 30 GHz t he ra dio spectrum flattens to a a 
0.65 (Mezger et al. 1986[ ) which possibly extends to the 
IR region if the flux of continuum emission measured 
at 6/xm by Tuffs et al. (1997) has a synchrotron origin. 

— Recent observations reveal hard X-ray emission that 
extends with a power-law photon index a x « —3 up to 



120 keV (Allen et al. |1997[ Favata et al. |1997| ); a syn- 
chrotron origin of this radiation implies acceleration of 
electrons to multi-TeV energies. 

Detailed studies based on high resolution mapping of 
Cas A at GHz frequencies have suggested that the iden- 



tification of bright compact radio features with electron 
acceleration sites can be problematic (ARLPB; Anderson 
& Rudnick 1996 , hereafter AR96). A particularly inter- 
esting finding of ARLPB and AR96 is a significant and 
rather unexpected correlation between the spectral index 
of the knots and their projected position in the shell, as 
well as their radio brightness: the steeper radio knots re- 
side mostly in the outer regions of the shell, at r ~ 2.5 pc, 
and tend to be brighter. 

The analysis in AR96 has shown that a self-consistent 
explanation of these correlations would be problematic 
if one assumes that the radio knots are the sites of ef- 
ficient acceleration of the electrons. For an SNR as young 
as Cas A the radiative energy losses cannot modify the 
spectrum of radio electrons. Therefore in the framework 
of a 'standard' spatially homogeneous source model ap- 
proach, one has to attribute the radio spectral indices ob- 
served to the source spectra of the electrons. The observed 
brightness/steepness trend of the radio knots would then 
apparently rule out effective electron acceleration in the 
knots, because such a process implies a hardening (see e.g. 
Berezhko & Volk 1997), rather than a steepening of the 
particle spectra, typically to the power-law index (3 ~ 2. 

There does exist however a natural way to modify the 
spectra of radio electrons, if we abandon the standard ap- 
proach of a spatially uniform source for treatment of these 
electrons. Spectral modifications become unavoidable if we 
take into account energy- dependent propagation and es- 
cape of relativistic particles from the regions of higher con- 
centration in an inhomogeneous medium. The efficiency of 
this process depends on the spatial gradients in the energy 
distribution N(r,E,t) of particles, and timescales of the 
spectral modifications can be as short as the escape time 
t(E). This effect is widely used for the interpretation of 
the galactic CR spectra in the framework of diffusive or 
Leaky Box models, but it has not yet been given proper 
attention in studies of radio emission in the CR sources 
themselves. 

In this paper we consider the consequensies of energy 
dependent propagation of relativistic electrons for the in- 
terpretation of the observed spectral, spatial and temporal 
characteristics of the broad-band nonthermal emission of 
Cas A. In Sect. 2 we introduce the 'two-zone' model for 
a spatially non-uniform radio source, separating compact 
regions with a high density of relativistic electrons (zone 
1) from the rest of the shell where the electron density 
is significantly lower (zone 2). In Sect. 2 we consider in a 
qualitative way possible consequensies of such an approach 
for the interpretation of the observed radio data. In Sect. 
3 we derive the system of kinetic equations for the electron 
energy distributions in the two zones. In Sect. 4 we assume 
that zone 1 components correspond to the sites of efficient 
electron acceleration, and show that this scenario is able 
to explain the broad band non-thermal radiation data of 
Cas A. In Sect. 5 we study the opposite scenario, which 
assumes that the enhancement of the electron density in 
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Fig. 1. VLA image of Cas A at A 6.3 cm, showing examples of morphological classes of radio features referred to in 
the text. Axes (tangential projection in RA, DE C for epoch B1950) are annotated at intervals of 1 pc (60.7 arcsec 
for the distance of 3.4 kpc derived by Reed et al. 1995|) , with the central (longer) tick indicating the B1950 ce ntre o f 
expansion of the system of fast optical filaments at RA 23:21:12.0 DEC +58:32:17.9 (Kamper & van den Bcrgh 1976 ). 
The greyscale spans the full dynamic range of 0.0 to 12.09 mJy per 0.375 arcsec pixel in 255 levels; a mild logarithmic 
scaling (0.85 exponent) has been applied to emphasise faint structures. The image was constructed from an almost 
fully sampled synthetic aperture using data from the VLA D, B, and A configurations at epochs 1983.0, 1984.0 and 
1983.8, respectively (see table 1 of Braun et al. (1987) for observational parameters). The brightness distribution has 
been deconvolved from the synthesised point spread function using the maximum entropy method (AIPS task VM). 
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zone 1 is caused not by active acceleration of particles 
there, but rather only by compression of the background 
population of relativistic electrons in the shell. We show 
that interpretation of the data within this latter scenario 
is problematic. In Sect. 6 we summarize the observational 
features which the model can explain, and in Sect. 7 we 
discuss implications and predictions of the model. 

2. The model 

2.1. Insufficiency of the single- zone approximation 

To lowest approximation Cas A might be represented as a 
homogeneous shell containing magnetic field, relativistic 
electrons and gas. Although not reflecting the real pat- 
tern of the source, this 'single-zone' approximation al- 
lows estimates of basic parameters such as the mean mag- 
netic field Bo and energetics in r adio e lectrons. In this 
approximation Cowsik & Sarkar (198C) have derived a 
lower limit to B > 8 x 1CT 5 G, comparing the expected 
bremsstrahlung flux of GeV electrons with the upper limit 
/( > 100 MeV) < 1.1 x 10- 6 ph/cm 2 s of SAS-2 (Fichtel et 
al. |1975Q and COS B detectors. 

With the recent upper limits of the EGRET tele- 
scope, /(> 100 MeV) < 1.2 x 10~ 7 ph/cm 2 s (Esposito et al. 
1996| ), this constraint on B is significantly strengthened. 



In Fig. 2 we present the fluxes of the bremsstrahlung and 
synchrotron radiations calculated for 3 values of the mean 
magnetic field: B Q = 10~ 4 G, 3.5 x 10~ 4 G and 7 x 10~ 4 G. 
For the mean gas density (in terms of 'H-atoms') we take 
riH = 15 cm~ 3 , corresponding to M ~ 15 Mq of matter in 
the shell (Fabian et al. 1980, Jansen et al. 1988| , Reed et al. 
19951 ), and 

use a mean atomic Z(Z + 1)/A = 4.3 derived 



by Cowsik & Sarkar ( [1980 ) from the elemental abundance 
estimates of Chevalier & Kirshner ( 1978| ) in Cas A. Con- 
tinuous injection of electrons into the shell, with a spec- 
trum Q(E) oc E~P° exp(-E/E c ) where O = 2a Q + 1 = 
2.54 and an exponential cutoff energy E c = 100 TeV, is 
assumed. The chosen E c explains the fluxes of hard X- 
rays above lOkeV as synchrotron emission for the case of 
Bo = 10~ 4 G. However, this value of Bo is unacceptable 
because the bremsstrahlung flux produced by GeV radio 
electrons is too high. The EGRET upper limit for the 7- 
ray flux of Cas A requires at least Bo — 3.5 x 10~ 4 G. Si- 
multaneously we have to assume a significant flattening in 
the source spectrum of electrons below 100 MeV in order 
to avoid contradictions with the observed X-ray fluxes at 
100 keV. For a single power-law distribution of electrons 
extending with fio — 2.54 down to the MeV region, the 
lower limit to the magnetic field is Bo = 7 x 10~ 4 G. Note 
however, that for Bo > 3.5 x 10~ 4 G the radiative energy 
losses of the electrons establish a spectral break at optical 
or lower frequencies, in which case the spatially homoge- 
neous single-zone model fails to explain the observed hard 
X-ray fluxes by synchrotron radiation. 
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Fig. 2. The fluxes of the synchrotron (heavy lines) and 
bremsstrahlung (thin lines) photons calculated in the 
framework of a spatially homogeneous single-zone model 
for the source function of accelerated electrons with power 
law index /3 = 2.54 for 3 different values of the mag- 
netic field: B = 10~ 4 G (solid), 3.5 x 10~ 4 G (dashed), 
and 7 x 10~ 4 G (dot-dashed). The radio flux measure- 
ments shown correspond to the data from Baars et al. 
(1977) corrected for the secular decrease to the epoch 
1986. The data at 1.3 mm (circle) and 1.2 mm (triangle) 
are from Chini et al. ( |l984f) and Mezger et al. ( |l986|) , re- 



spectively. The hatched region correspo nds t o the range 
of X-ray fluxes measured by Allen et al. ( 1997 ). Note that 
the RXTE fluxes are shown only above lOkeV since below 
lOkeV the X-ray fluxes are dominated b y the rmal emis- 
sion. The OSSE measurements (The et al. 1996 ) are shown 
by crosses, and the bar corresponds t o the EGRET upper 
limit above 100 MeV (Esposito et al. [1996D. 



2.2. The two- zone model 

A more realistic model for Cas A has to account for the 
inhomogeneities in the radio brightness distribution and 
for the observed spread in the power-law exponents of the 
radio spectra of individual structures within a ~ (0.6 — 
0.9). 

The simplest approximation for a non-uniform radio 
source corresponds to a 'two-zone' model, where the source 
with total volume Vo consists of 2 different regions, or 
zones, with volumes Vi and V 2 = Vo — V\, with 2 dif- 
ferent spatial densities (concentrations) of radio electrons 
ni {E) = N X {E)/Vi and n 2 {E) = N 2 {E)/V 2 . Propagation 
effects may result in a significant difference between the 
two-zone and single-zone models only if these densities 
are essentially different. Therefore the principal definition 
for these two zones is that the concentration of electrons 
in zone 1 is much higher than in zone 2, ni 3> n 2 . The 
energy distributions of the electrons in these zones can 
be different. Therefore in the power-law approximation 
N\, 2 {E) oc E -131 - 2 the exponents (3\ and f3 2 can be dif- 
ferent. Because the spectral ratio n\{E)/n 2 (E) may thus 
depend on energy, it is convenient to specify that the con- 
dition ni ^> n 2 relates first of all to electrons with energies 
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E r*> 1 GeV which are typically responsible for radio emis- 
sion in the '1 GHz' band. 

In order to apply the two-zone model to Cas A, we have 
to foresee some observational manifestations of the radio 
emission of zones 1 and 2. It is reasonable to expect that 
not only the electron densities, but also the mean mag- 
netic fields B\ and in the two zones differ, with the 
magnetic field presumably being higher in the regions of 
higher electron densities. Therefore, the radio emissivities 
in these zones will be strongly contrasting, with regions 
belonging to zone 1 characterized by a much larger emis- 
sivity than regions in zone 2. It is also plausible that zone 
1 components should be rather compact, with a small vol- 
ume filling factor in the source, or otherwise the overall 
emission would be strongly dominated by this single zone. 

Application of this criterion to Cas A immediately 
places all bright fragments of the radio ring in zone 1. 
Other constituents of zone 1 are the radio knots (and bow 
shocks) most of which are characterized by a strongly en- 
hanced cmissivity. We note however that in reality the 
knots embrace a wide range of emissivities, and because 
the simplified two-zone approach allows only 2 different 
electron densities (and emissivities), some compact struc- 
tures with relatively low brightness should be more cor- 
rectly classified as belonging to zone 2. This may have a 
physical basis - enhanced emissivity due to mild compres- 
sion of magnetic fields, but without strong excess in the 
electon densities. Zone 2 then consists of the rest of the 
(spherical) shell between the radio ring and the outer edge 
of Cas A, i.e. the main body of the radio plateau, but also 
includes some low-brightness radio knots. 

For purposes of practical classification of low-intensity 
knots it may prove useful to resort also to the second po- 
tentially observable distinction between zones 1 and 2 - 
different mean spectral indices a\ and ui of their intrinsic 
emission. Due to a faster escape of higher energy electrons 
from the regions of high concentration, one could gener- 
ally expect that the radiation spectra produced in zone 1 
should be steeper, ol\ > (see below Sect. 2.2.3). Because 
the spectral index of the total emission a ss 0.77 is to be 
maintained, an immediate conclusion is that the radiation 
spectrum of the diffuse zone 2 should be noticeably flatter, 
than the mean, i.e. cti < a® < ct\. 

Because of the complicated morphology, a quantita- 
tive treatment of the spectral indices of different types of 
structures is problematic in practice. In order to assess 
the characteristic values for 0:1,2, we note that the line of 
sight integrated radio brightness seen towards individual 
radio knots shows a broad range of spectral indices from 
a m i n ~ 0.6 to a max ~ 0.95, as measured by AR96 be- 
tween 1.4 GHz and 5 GHz. These spectral indices however 
do not always correspond directly to the intrinsic spectral 
indices of the knots, due to background contamination 
from the plateau. Although the detailed least square fit- 
ting procedure for spectral index measurements used in 
AR96 (their Sect. 2.4) largely removes this problem for 



knots with high contrast to the plateau, for faint knots 
at the limit of confusion with the plateau, the spectral 
index will be biased towards the (presumably flatter) in- 
dex of that emission. In this regard, the result of AR96 
presented in their Fig.6d appears as very informative. It 
shows that most of the radio knots found at large angular 
distances 8 > 130" (the projected r± = 2.17pc) belong 
to the steep-spectrum population: from 62 such knots in 
total, 45 have indices a > 0.8, and 22 (i.e. 35 % !) be- 
long to the 'extreme' end of the observed spectral indices 
0.86 < a < 0.95. Because at these projected distances 
the brightness of the plateau is generally reduced]], this 
may indicate that the intrinsic spectral index of the radio 
knots on average could be much steeper than the mean 
ao, and that the spectral index a > 0.86 of the 'extremely 
steep' population of the observed knots might be not so 
extreme, but rather a representative value for the intrinsic 
index of the knots. An (indirect) indication of this sug- 
gestion could be also another result shown in Fig. 6d of 
AR96: 83 from 123 knots found at angular distances close 
to the radio ring, 80" < 8 < 120", where the brightness of 
the plateau is relatively high, belong to the population of 
knots with an 'intermediate' steepness 0.71 < a b s < 0.8. 
As discussed below in Sect. 6, the model does not exclude 
a physical reason for such geometrical correlation between 
the knot index and its proximity to the ring. However, this 
effect can be explained also as a result of a line-of-sight 
mixture of comparable fluxes of the steep-spectrum knots 
with intrinsic a ~ (ao + 0.1) and of the flat-spectrum dif- 
fuse plateau emission with a ~ (0.6 — 0.65) (see this Sect, 
below). Moreover, if we believe that the real geometry of 
the shell is quasi-spherical, it is difficult to suggest any 
other reasonable explanation, except for invoking the flux 
contamination effect, why only 2 knots from those 123 
show the index a bs > 0.83, whereas the expected num- 
ber of 'extremely steep' knots to be located at deprojected 
radii r > 2.17pc but at the angular distances close to the 
"ring" should have been comparable with ^ the number 
22. 

We can conventionally define that all knots which have 
a moderate brightness and show spectral index a < 0.7 be- 
long to zone 2. Then, given the tendency for brighter knots 
to be steeper (AR96), this classification would retain most 



We note that the low frequency radio maps of KPDE, with 
angular resolution too poor to see individual radio knots, show 
a profound gradient of steepening spectral index of the plateau 
towards a < 0.9 with projected radius r > -Rring- It is plau- 
sible (though not quantitatively investigated) that this trend 
is due to the population of steep spectra radio knots detected 
by AR96, which at lower angular resolution would contribute 
more flux than the plateau at angular distances > 130"; the 
intrinsic spectrum of the plateau emission might be then flat. 

2 Even if the steep spectrum knots were physically dis- 
tributed on the outermost radius of Cas A (~ 150"), we would 
still expect significantly more than 2 knots to be seen within 
projected distances 80" < 6 < 120" 
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of the radio knots (at least in terms of their overal flux) in 
zone 1. At the same time, excluding thus the Hat-spectrum 
knots from zone 1, for the average intrinsic spectral index 
of the remaining population of knots one can reasonably 
suppose a^tn ~ (a + 0.1) = 0.87 or so. 

It is possible that the intrinsic spectral index of the 
zone 1 components identified with the bright fragments of 
the radio ring is also close to this value. Th e low frequency 
measurements of Woan & Duffet-Smith (1 



In order for the mean spectral index between = 1.4 GHz 
and z/** = 5 GHz to be equal to ao, one needs a flux ratio 



199C: hereafter 



WDS90) at angular resolution > 14" show that the av- 
erage spectral index of the bright large radio components 
correlated with the ring is ~ 0.83, i.e. again significantly 
larger than the mean ao of the total emission. If we accept 
that the contribution of the diffuse emission, zone 2, may 
be flat, the intrinsic spectral index of the fragmented ra- 
dio ring should be even larger. Then ai ~ (ao + 0.1) may 
be a reasonable value for the average spectral index of the 
entire zone 1. For qualitative discussions in this section we 
will use ot\ — 0.87. We note however that our model does 
not pretend to predict characteristic spectral indices with 
the accuracy of the second digit after the comma. 

The first zone will thus include a number K ~ 200 in- 
dividual components, i.e. basically the bright fragmented 
radio ring, the 'radio' bow shocks, and the main part of 
all radio knots. Besides, zone 1 could also contain some 
regions of high concentration of electrons and high emis- 
sivities which are not immediately distinguished on the 
intensity maps either because of their very small sizes or 
because of strong contamination by, and confusion with 
the diffuse plateau emission. The latter may be the case 
if we believe that the fragmented "radio ring" actually 
represents a chain of large but geometrically flat (thin) 
structures of the "radio sphere" which are observed mostly 
'edge-on'. Then at angular distances interior to the ring, 
9 < 80" or so, such fragments observed 'face-on' may eas- 
ily merge with the plateau emission and show up as diffuse 
'clouds' of only moderately enhanced brightness (see Fig. 
I)- 

The flux produced in zone 1 can be estimated if we 
remember that at 5 GHz the (background subtracted) flux 
of the radio ring fragments makes up w 30%, and that 
radio knots produce w 20% of the total flux (Tuffs |1986| ). 
The latter value should be somewhat reduced for the flux 
of knots to be removed from zone 1 and classified in zone 2. 
On the other hand, it seems quite possible that up to (5 — 
10) % of the total flux should be assigned to the 'unseen' 
population of zone 1 (especially, to face-on structures of 
the "radio sphere"). Therefore, we can estimate the overall 
flux Ji of zone 1 as ~ (45 — 55) % of the total, concluding 
that the fluxes of zone 1 and zone 2 at 5 GHz are about 
the same, Ji ~ J 2 . 

The spectral index a 2 in the diffuse zone 2 can be then 
predicted if we note that the total flux has the form 



,h _ (iW») ai ~ a ° ~ 1 

Ji ~~ i - (iW*) a2_ao 



(2) 



J(u) = Ji 



,5 GHz 



5GH2 



(1) 



Then for J 2 — J\, olq ~ 0.77, and a.\ — 0.87, we find 
a 2 = 0.65. It may somewhat vary assuming slightly differ- 
ent J 2 1 J\ or ol\ (or ao , which can be also slightly less than 
0.77, see Fig. 2). In particular, a 2 may be rather close to 0.6 
if we allow for the idea that the 'unseen' face-on structures 
of the "radio sphere" could contribute into zone 1 a third 
of what is contributed by the edge-on "radio ring". Re- 
markably, these values of a 2 correlate well with the lower 
range of spectral indices of the flat-spectrum radio knots. 
This suggests that flat-spectrum radio knots, with spec- 
tral indices close to 0.6 (or flatter than 0.7 as classified 
above) may indeed represent the regions of moderate lo- 
cal compression of the magnetic field and/or relativistic 
electrons of the diffuse shell, and therefore may give infor- 
mation about the spectrum of the background population 
of electrons there. 

Thus, the model suggests a rather flat intrinsic spec- 
trum of the radio plateau in the range a 2 — 0.6 — 0.65, 
which corresponds to the spectral index of electrons f3 2 ~ 
2.2 — 2.3. For zone 1, a\ ~ 0.87 implies a rather steep elec- 
tron spectrum with /3i ~ 2.74. Taking into account the 
enhanced brightness of the compact zone 1, we will refer 
to zone 1 structures as the compact bright steep-spectrum 
radio (CBSR) components. 

2.2.1. Synchrotron self- absorption in the two-zone model 

Even though still highly simplified, this two-zone approach 
suggests a qualitative interpretation of several observed 
radio features of Cas A. In particular, it allows us to solve 
the problem of the low- frequency turnover below 20 MHz. 

Although the measurements of KPDE at ~ 100 MHz 
show evidence for thermally absorbing low-temperature 
gas in the central < 1 pc region of Cas A, the thermal ab- 
sorption cannot result in any noticeable effect in the shell 
where T ~ 10 7 K. On the other hand, in the case of a spa- 
tially homogeneous source (the shell) with a mean field 
Bo = 0.7 mG, the process of synchrotron self-absorption 
results in a turnover of the radio spectrum only at fre- 
quencies below 5 MHz (see Fig. 2). To explain the observed 
turnover below 20 MHz, in a single-zone approach one 
needs Bq > 5 mG. But this assumption boosts the mag- 
netic field energy up to Wb > 10 51 erg, and therefore is 
hardly acceptable. 

The situation essentially changes if one takes spatial 
inhomogeneities in the distribution of radio electrons and 
magnetic fields into account. Indeed, in the framework 
of the proposed two-zone approach the radio flux at fre- 
quencies well below 1 GHz is dominated by the steep- 
spectrum component, i.e. one has to require synchrotron 
self- absorption mainly of the zone 1 flux at v < 20 MHz. 
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Because of the very high density of the radio electrons and 
of the magnetic field in compact zone 1 components, syn- 
chrotron self-absorption there should be expected at fre- 
quencies significantly higher than in the single-zone model. 

In order to estimate the model parameters needed 
for synchrotron self-absorption to occur below 20 MHz, 
we note that for a power-law distribution of electrons 
N(E) oc E~P the synchrotron luminosity can be expressed 
as 



L{v) = 2 x 1CT 25 Cp N* 



D 



lmG 



1+0 



10 GHz 



1-/3 



erg 
sHz 



(3) 



where N* = E^N{E, f ) with = 1 GeV. Cp is a weak 
function of P that changes within 1 ±0.15 for /? ~ (2 — 3), 
and we will assume Cp = 1. With the same accuracy, the 
synchrotron absorption coefficient is given by 



k{v) ~ 9x 10" 17 /310 3 ^ 2 ^ / B 



V VlmG 



1+0/2 



V \-(2+/3/2) 



(l0 7 Hz) 



(4) 



where V is the volume in units of cm . The volume of 
zone 1 can be represented in the form V\ = Y2iLi — 
S pl d, where d is the mean thickness of individual CBSR 
components along the line of sight and S'pr = Sj=i 
is their overall surface in the plane of the sky. Requiring 
the synchrotron opacity to be significant at v < 20 MHz 
(which implies nd ~ 1 at v ~ 15 MHz) we obtain 



S m ~ 1.3 x 10 37 ' J[ 



B 1 



400 Jy/ VlmG 



0.5 



(5) 



For the external radius of the plateau region R = 2.5 pc 
this leads to an area filling factor 



So 



0.07 



Ji 



400 Jy/ VlmG 



5, 



0.5 



(6) 



where So = ttRq- For B\ ~ (1 — 2) mG in zone 1 (see Sect. 
4), the ratio S pr /So ~ 0.1. From Fig.l it is apparent that 
the area filling factor of all bright radio features is indeed 
of the order of 10 % . 

From Eq.(5) it follows that the characteristic length 
scale for K ~ 200 CBSR components of the zone 1 should 
be of order Al ~ (S pT / ' K) 1 / 2 ~ 0.1 pc. If we separate the 
radio knots from the larger fragments of the radio ring, 
and take into account that the flux J^ n of the knots con- 
stitutes about a third of J\ , the mean radius of the knots 
can be estimated as ~ 0.03 pc. At the distance 3.4 kpc this 
corresponds to the angular sizes of about 2", very similar 
to what is really observed. 



2.2.2. Flattening of the radio spectrum 

The superposition of flat and steep spectral components, 
with Ji ~ J2 around 5 GHz, results in a flattening of the 
total radio spe ctrum at higher frequencies. Observations of 
Mezger et al. ( 1986 ) do indicate a noticeable flattening of 
the Cas A spectrum in the wavelength range from 1 cm to 
1mm, with power law index ~0.65 that possibly extends 
up to the near-IR region if the flux of continuum emission 
measured at 6 /im by Tuffs et al. (1997) is synchrotron in 
origin (see Fig. 2). 

Another important effect related to such composite 
fluxes consists in the possibility to give a new interpre- 
tation of the claimed secular flattening of the radio spec- 
tra (Dent et al. [1974^ O'Sullivan & Green |199S| ). Namely, 
the flattening of the radio spectra in time could be easily 
explained as being due to different rates of decline of the 
fluxes of the two zones, with the plateau emission drop- 
ping more slowly. Observations at 5 GHz have shown that 
the CBSR components are fading more rapidly than the 
total emission, which then requires that the plateau must 



be fading less rapidly (see Table 6.1 of Tuffs 1983). If the 



plateau indeed were to have a flatter spectral index than 
the CBSR components, this would result in a flattening 
of the spectral index of Cas A with time. Note that pre- 
viously, in the framework of the single-zone approach, the 
secular flattening could be explained only by the assump- 
tion of a gradual flattening of the spectrum of electrons 
eithe r at higher ene rgies o r in time (Scott fc C hevalier 

" Cowsik & Sarkar 11984 Ellison 



1978 



1975 , Cheva lier et al. 
& Reynolds |l99l| ). 



2.2.3. Effects of energy dependent propagation 

Another advantage of a spatially inhomogeneous model is 
its ability to explain the basic characteristics of the ob- 
served nonthermal radiation in terms of very efficient ac- 
celeration processes operating in Cas A, which generally 
result in a source function of accelerated particles with 
hard spectral indices /3 acc ~ (2 — 2.2) (e.g. efficient Fermi 
acceleration at strong shock fronts) . Even though the spec- 
tral index formally corresponding to the mean ao — 0.77 
is /?o > 2.5, the assumption of /3 acc much harder than 
Po becomes possible if one takes into account the energy 
dependent propagation and escape of electrons from re- 
gions of high concentration on timescales T esc (E) oc E~ s . 
For example, assuming that the sites of efficient parti- 
cle acceleration are located in zone 1, one could expect 
that the electron energy distribution N\{E) formed there 
on timescales t t csc would have the power-law index 
Pi = Puce + S. On the other hand the leakage from zone 1 
corresponds to injection of relativistic electrons into sur- 
rounding zone 2 with the rate Q%{E) — Ni/t csc oc E~P' 2 
where P% = P\ — 8 = /3 acc . With a reasonable source func- 
tion index /3 acc ~ 2.2 — 2.3 and energy-dependent escape 
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with S ~ 0.5, the values of /3i and fa needed for interpre- 
tation of the Cas A radio data could be reproduced. 

Thus, the energy distribution of electrons Ni (E) in the 
compact regions of efficient acceleration, which results in 
the enhanced density of relativistic electrons, may be actu- 
ally much steeper than N 2 (E) in the surrounding volume. 
Interestingly although the electron distribution N 2 (E) in 
the plateau region would show the hard spectral index of 
acceleration /3 acc , in principle no acceleration of electrons 
needs occur there at all. It is important to note, however, 
that on a qualitative level the same power-law indices 
02 and /?i for both zones could be expected also in the 
case when relativistic electrons were accelerated in zone 

2, while the electron density in CBSR structures would be 
enhanced due to local compression of the electrons, but 
not efficient acceleration in zone 1. To distinguish between 
these 2 basic scenarios, quantitative modelling is needed. 

3. Kinetic equations in the two-zone model 

Approximating the momentum distribution function of 
relativistic electrons as an isotropic function f p of the 
momentum p = p | , the kinetic equation for the energy 
distribution / = f(r,E,t)dE — Airp 2 f p dp of relativistic 
particles (E ~ pc) can be written as: 

d -l = div r p g rad r /)-div r (u/) + ^A (jB/) 



(7) 



Here u = u(r, t) is the fluid velocity, P = P(r, E, t) is the 
energy loss rate of electrons. Assuming isotropic diffusion, 
D = D(r,E,t) is the scalar spatial diffusion coefficient, 
whereas A[f] is a functional standing for various stochastic 
acceleration terms of electrons. 

In the two-zone model the source with volume Vo is 
subdivided into zone 1 composed of K compact structures, 
with volume V\ = Y^f=i ^ Vo, spatially immersed in 
the zone 2 with Vz = Vo — Vi « Vo ■ The set of equations de- 
scribing the evolution of the total energy distribution func- 
tion of particles in these two zones, Nj(E,t) = J f dVj , 
can be found by integration of Eq.(7) over Vj (j — 1,2). 
Integration of the term df jdt over V\ results in dNi/dt. 
The volume integral of the two first terms in the right hand 
side of Eq.(7) is reduced to the sum of K integrals on the 
surface S\2 = "Y^f Si separating the zone 1 components 
from zone 2: 



(8) 



/ [div r (Pgrad r /)-div r (u/)]d 3 r = 

V [ J> £>(egrad r /)ds - <f> (eu)/ds] . 

~[ JS Z JSi 



Here e is the unit vector perpendicular to the surface el- 
ement ds, directed outward from zone 1 to zone 2. These 



terms describe the rate of diffusive and convective ex- 
change of particles between the two zones. 

To proceed further, let us separate those parts of the 
interface S\ 2 where the scalar product (u e) has positive 
and negative signs: Si 2 = S + + S-. The surface S + corre- 
sponds to the regions of plasma outflow from zone 1 into 
zone 2, so the integral over S + of the second term in Eq.(8) 
describes the convective escape of electrons, —Ni(E,t)/T c , 
where 



h/ui , 



(9) 



h is characteristic thickness of, and U\ is the plasma ve- 
locity in CBSR components. 

The first integral in Eq. (8) over the surface S+ can be 
simplified if we approximate (egrad/) ~ (f 2 — fi)/Al, 
where /1 = f\(E,t) and f 2 — f 2 (E,t) are the volume 
averaged distribution functions of relativistic particles in 
zones 1 and 2, respectively, and Al = Al(E, t) is the char- 
acteristic thickness of the transition layer between these 
zones. Taking into account that the volume of a body can 
be expressed as V — ahS + , where h is the thickness and 
a ~ 1 is a factor depending on the shape of the body, we 
find 



K K 

W £(egrad/) = ]T 

1=1 JSi+ i = l 



h - h 



Al 

h - h 



Di Si-\ 



Tdif 



V 



(10) 



where Tdif has the meaning of characteristic diffusive es- 
cape time of relativistic electrons from the regions of 
higher electron densities: 



1 = y> Di Vi (t) ) D T (E,t) 
T dii (E,t) j^ahiAk Vi ~* ahAl(E,t) 



(11) 



Here D T (E,t) corresponds to the mean of the spatial dif- 
fusion coefficient D r (r, E, t) on the surface S+. Taking fur- 
ther into account that Nj = Vj fj (j = 1,2), the diffusive 
escape (exchange) term is reduced to the form 



ViN 2 (E,t) N^E,t) 
V 2 T di {(E,t) r dil (E,t) 



(12) 



The integrals over the surface 5_ in Eq.(8) describe 
the inflow of electrons from zone 2, with some rate Q21, 
through the frontal surfaces of zone 1 components. For a 
mainly adiabatic compression of relativistic particles (see 
Sect. 5), Q21 could be directly connected with the energy 
distribution N 2 of electrons in zone 2. One could how- 
ever expect the formation of strong shocks due to fast mo- 
tion of CBSR components with respect to the surrounding 
plasma, so that the upstream electrons could be acceler- 
ated before entering zone 1. In that case Q21 would cor- 
respond to the source function of electrons accelerated on 
the shock fronts and injected into zone 1. 
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An essential part of the diffusive shock acceleration 
process is given by the adiabatic compression term in 
Eq.(7) (third term) in the shock region near the surface 
S- . Therefore the volume integral of that term in a thin 
region adjacent to S- is to be formally incorporated into 
the source function Q2 1 . The integral in the main part of 
the volume V\ describes adiabatic energy losses (or gain) 
of electrons, and can be combined with the volume integral 
of the forth term in Eq.(7) as: 



d_ 
~dE 



(£div r u/3 + P) fd 3 r 



V x 



d_ 
~dE 



(13) 



where Pi = Pi(E,t) is the mean rate of energy losses of 
electrons, which now includes also the adiabatic energy 
loss term P ad ~ E $(ue) ds/3. 

Finally, the volume integral of the last term in Eq. (7) 
describes internal sources of accelerated electrons Q^ nt ^ 
inside the volume V\. The final equation then reads: 



dt 



dE 



V 2 T di f 



(14) 



Here Qi = Qi(E,t) = Q 2 i + Qi corresponds to the 
total injection rate of electrons in the first zone, and r csc 
is the overall " diffusive + convective" escape time of par- 
ticles 



1 



f 



Tdif CM T&) 



(15) 



Because Tdif generally increases for decreasing E, T esc 
becomes energy independent at sufficiently small energies 
when Tdif > t c . The range of actual energy dependence of 
t csc is limited also at high energies since the diffusive es- 
cape time cannot be less than the light travel time across 
the CBSR components ~ h/c. This obvious requirement 
formally follows from the condition that the characteris- 
tic lengthscale Al(E) for spatial gradients in the distri- 
bution function f(r,E,t) cannot be less than the mean 
electron scattering path \ SC (E), since otherwise the dif- 
fusion approximation implied in Eq.(7) fails. Taking into 
account that D ~ A sc c/3, from Eq.(ll) we find that in- 
deed Tdif(-E,t) > T m i n ~ 3ah/c. Assuming D(E) oc E Sl , a 
reasonable approximation for the diffusive escape time is 
then 



r dif (£,t) =r,(t) (E/E.) - s + T min (t) 



(16) 



From Eq.(ll) follows that generally the power-law in- 
dex 8 for diffusive escape time should be smaller than the 
index Si of the diffusion coefficient: since faster diffusion of 
more energetic particles tends to smooth out the gradients 
of the distribution function /(r, E, t) more effectively, the 
characteristic lengthscale Al(E, t) should increase with en- 
ergy. Then for a power-law approximation AZ oc E S2 the 
index 62 > 0, so 8 = 61—62 < 8%. This conclusion is 
important since it allows us to assume a power-law index 



6 ~ 0.5 — 0.6 needed for modification of the electron spec- 
tral index, not excluding at the same time efficient shock 
acceleration of particles in the Bohm diffusion limit which 
corresponds to #1 = 1. 

Integration of Eq.(7) over the volume of zone 2 results 
in an equation for N2{E,t): 



8N 2 d( P 2 N 2 ) 



Ot 



dE 



^CSC 



vlN 2 

V 2 T di f 



Q-2 



(17) 



where P2 is the mean energy loss rate, and Q2 is the source 
function of electrons in zone 2. The second and third terms 
here result from the same surface integrals on the inter- 
face between zones 1 and 2 as in Eq.(8), but they have 
opposite signs since now the direction of the unit vector 
e is reversed. Note that formally the volume integrals in 
Eq.(8) result in two additional surface integrals describing 
particle exchange between zone 2 and the exterior on its 
outer boundary. However, this actually extends the two- 
zone model to a three-zone one, therefore we neglect here 
these exchange terms in order to remain in the framework 
of the two-zone approach. 

The set of kinetic equations (14) and (17) describes 
the evolution of relativistic electrons in the general case 
of an expanding inhomogeneous source when the electron 
energy losses and particle exchange rates between two 
zones are time-dependent. Given these parameters and 
the source functions Qi^,{E,t), the electron energy dis- 
tributions can be found numerically by iterations, using 
analytic solutions to the kinetic equation of electrons in 
an expanding homogeneous medium (Atoyan & Aharonian 
1999 ). To reduce the number of free parameters, we shall 
consider basically the case of a stationary source because 
the current energy distributions depend mostly on the in- 
jection of the electrons during recent (largest) timescales, 
A^o — 300 yr. Only for calculations of the secular decrease 
of the fluxes we shall use the general results for an ex- 
panding source, taking into account adiabatic losses and 
decreasing average magnetic fields. 

4. CBSR components as electron acceleration 
sites 

Let us first consider the hypothesis that the compact 
bright radio components are the main sites of electron 
acceleration, and neglect possible acceleration in zone 2, 
i.e. Q2 — 0. The electron injection function in zone 1 is 
approximated as 



Q(E) = Qr,E-^exp(-E/E, 



(18) 



where we have omitted the suffix T' in Q. For calculations 
we assume equal radii Ri and magnetic fields Bi for zone 
1 components. 

Since the thickness ft of a spherical object corre- 
sponds to its diameter, Eq. (9) gives t c = 2Ri/m. For 
Ri ~ (0.05 — 0.1) pc, and u\ less than the sound speed 
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v s — ■J ^adbkT/Amp ~ 150km/s in the gas with T\ < 
3 x 10 7 K and mean atomic A ~ 15 (corr esponding to 
Z(Z+l)/A ~ 4.3, Cowsik & Sarkar |1980| ), the time r c 
is larger than the source age to ~ 300 yr. Therefore the 
precise value of u\ will not significantly affect the results. 
The minimum escape time of electrons from the zone 1 is 
taken as T m - m = R\/c. 




Fig. 3. The energy distributions of the electrons in zone 
1 (solid) and zone 2 (dashed), calculated in the two-zone 
model assuming that the CBSR components in CAS A 
correspond to the sites of electron acceleration. The dot- 
dashed curve corresponds to the initial iteration for the 
energy distribution of electrons in zone 1 when the 're- 
turn' diffusive flux of electrons from zone 2 into zone 1 is 
neglected (the term oc A^/Tdif in Eq.14). Parameters for 
the source function are /3 acc = 2.24, E c = 3.5 x 10 13 eV; the 
magnetic fields B x = 1.5 x 10~ 3 G and B 2 = 4x 10~ 4 G; 
the diffusive escape time is defined by = 25 yr and 
S = 0.6. For zone 1, K = 200 components with the radius 
i?i = 0.05 pc are supposed. 

The time r* in Eq.(l) defines the ratio of the mag- 
netic fields in two zones needed for an explanation of the 
flux ratio J\/J% ~ 1. Indeed an effective modification of 
the injection spectrum of radio electrons in the GeV re- 
gion implies that at these energies the overall escape time 
Tesc(E) ~ t*(E/E*) s , and that it is significantly smaller 
than to ~ 300 yr. The leakage of the electrons from V\ 
will result in N2(E) ~ Ni(E) to / t csc (E) electrons accu- 
mulated in the surrounding zone 2. Using also Eq.(3), we 
then find 



Bi ( t 
Thus, for r, 



2B X 
lmG 



1 + 02 



(19) 



< 0.H o the ratio B 1 /B 2 > 3. 
In Fig. 3 we show the energy distributions of electrons 
ni,2(E) = Ni t 2{E,to)/Vi£ calculated the injection spec- 
trum with /3 acc = 2.24 and E c = 35TeV, and the escape 
time parameters t* = 25 yr and 6 = 0.6. At GeV ener- 
gies the power law index of electrons in zone 1 is indeed 
strongly modified, almost reaching the maximum theo- 
retical value Pi » (3 acc + S. However both at lower and 



higher energies there is some flattening of n±(E). At en- 
ergies E < 100 MeV the flattening is connected with the 
increase of the escape time to raif > 100 yr, which be- 
comes comparable to the age of Cas A, reducing thus the 
efficiency of the ecape. At high energies t^h (E) is limited 
by the energy-independent minimum escape time r m i n in 
Eq.(16). This effect can be seen in Fig. 3 for the dot-dashed 
curve at E > 100 GeV. Note, however, that for the self- 
consistently calculated spectrum (solid curve) of electrons 
in zone 1 the flattening is noticeable already at energies 
E > 10 GeV where n\(E) becomes comparable to n 2 {E). 
At these energies the gradients in the spatial distribution 
of relativistic electrons in the source become small, reduc- 
ing the net exchange of particles between the zones, which 
again results in a less efficient modification of Q(E). 
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Fig. 4. The fluxes of synchrotron (heavy lines) and 
bremsstrahlung (thin lines) photons produced by electrons 
shown in Fig. 3. The dashed and dot-dashed curves corre- 
spond to the fluxes expected from the compact radio struc- 
tures and the diffuse plateau, respectively, and the solid 
curves correspond to the synchrotron and bremsstrahlung 
fluxes of two zones. The full dots show the total nonther- 
mal flux. 
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Fig. 5. The spectral indices of the synchrotron fluxes cor- 
responding to the total (solid), zone 1 (dashed), and zone 
2 (dot-dashed) radiations shown in Fig. 4 
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The synchrotron and bremsstrahlung fluxes are pre- 
sented in Fig. 4, and Fig. 5 shows the spectral indices at 
different frequencies. Because the synchrotron radiation 
in zone 1 is steep, it becomes dominant at low frequen- 
cies. The compactness of zone 1 components leads to syn- 
chrotron self absorption at v < 20 MHz. At frequencies 
above 5 GHz the plateau emission dominates resulting in 
a noticeable flattening of the overall spectrum in the far 
infrared (FIR). All these results are in agreement with 
expectations derived qualitatively in Sect. 2. 

For the model parameters used in Fig. 4 the energy 
content in relativistic electrons in zones 1 and 2 is W e i = 
1.3 x 10 48 erg and W c2 = 1.4 x 10 48 erg, respectively, and 
the energy in the magnetic fields Wbi = 2.6 x 10 47 erg and 
Wb2 — 6.8 x 10 48 erg. The total energy in relativistic elec- 
trons, W e = Wci + Wo2, can be estimated analytically tak- 
ing into account that: (a) the shape of the overall energy 
distribution of electrons N(E) = N^E) + N 2 {E) repeats 
the power law injection oc E~P* CC until the break energy 
~ 300 GeV due to synchrotron losses (in zone 2), and (b) 
the total number of GeV electrons TV, = E* N(E^) (with 
E* = 1 GeV ~ 2000 m c c 2 ) is defined mainly by electrons 
in zone 2, « N*2, which can be estimated using Eq.(3). 
For /3 acc > 2.1 we find 



W, 



2.2 x 10 46 



(2000) 



/3a 



B 2 



erg. (20) 



/? acc -2 \lrnG, 

For /3 acc = 2.24 and B 2 = 0.4 mG this equation results 
in W e w 2.5 x 10 48 crg, in agreement with the numerical 
results. 




Fig. 6. Radiation fluxes calculated for /? acc = 2.2 and 3 
different magnetic fields B\ in zone 1, B\ = 1 mG (solid), 
1.6 mG (dashed), and 3mG (dot-dashed), assuming the 
same ratio B\jBi = 4. Exponential cutoff energies E c 
and escape times r* in these 3 cases are: E c = 16 TeV and 
r* = 26 yr (solid), E c = 27 TeV and t* = 20 yr (dashed), 
and E c = 200 TeV and r* = 14 yr (dot-dashed). For the 
assumed here number of zone 1 components K = 150 with 
radius R\ — 0.06 pc the sky covering fraction of the zone 
1 is 0.09. 

In Fig. 6 we show the fluxes calculated for /3 acc = 2.2 
and different values of the magnetic fields in zones 1 and 



2. For all 3 cases the spectra of synchrotron radiation co- 
incide up to the IR, but at higher frequencies they differ 
because of radiative break induced by synchrotron losses. 
Note that for interpretation of the X-ray fluxes above 
lOkeV one needs higher values of the exponential cut- 
off energy for higher magnetic field B\: E c = 16 TeV 
for B x = lmG and E c = 27 TeV for B 1 = 1.6 mG. 
For Bi = 3 mG the calculated synchrotron fluxes remain 
significantly below the observational data even assuming 
E c = 200 TeV. 

However, the synchrotron losses in high magnetic fields 
effectively limit the maximum energy of electrons to val- 
ues well below E c ~ 200 TeV. For example, in the case 
of diffusive shock acceleration, the acceleration rate can 
be estimated as dE/dt ~ Eu\/D{E), where D(E) is the 
diffusion coefficient and u 2 is the speed of CBSR compo- 
nents with respect to the surrounding plasma. Equating 
this rate in the Bohm diffusion limit to the synchrotron 
loss rate, the characteristic maximum energy of electrons 
can be estimated as: 



£(max) _ 3 x 1Q 



13 



U2 



3000 km/s/ VlmG 



B 



-1/2 



eV. (21) 



Comparison with the values of E c in Fig. 6 shows that 
diffusive shock acceleration could be consistent with syn- 
chrotron origin of the hard X-rays from Cas A for B\ < 
2mG. Higher values of B\ would require a more effi- 
cient acceleration mechanism. But anyway, the value of 
Bi = 3 mG is a very conservative upper limit of our model 
for the mean magnetic field in the bright radio structures 
if the measured hard X-ray flux has synchrotron origin. 
At the same time, the lower limit for B\ is about lmG, 
which is needed for consistency of the bremsstrahlung with 
the flux upper limits of EGRET and OSSE detectors (see 
Fig. 6). It is interesting that magnetic fields corresponding 
to the energy equipartition with relativistic electrons in 
the CBSR components are in the range B\ ~ (1.5 — 2) mG. 

Another feature to be addressed in this section is the 
wavelength-dependent rate — din J(v)/dt of the secular 
decline of radio fluxes. For calculations of this rate we 
need to specify the rates of the magnetic field decline and 
the adiabatic losses of electrons. For a spherical source 
(shell) homogeneously expanding with speed uq, the adi- 



abatic losses P a d = E/t ao \2 , where t ao \2 — t c 



Ro/u 



1000 yr corresponds to the current expansion time ('age') 
of the outer shell. In the 'standard' power-law approx- 
imation for the magnetic field declining with the shell 
radius as B(R ) oc Rq m , the characteristic decline time 
t-g, = B/(—dB/dt) = t cxp /m. For a magnetic field 'frozen' 
into the spherically expanding source m = 2, but m < 2 if 
the fields are effectively created in the shell. The assump- 
tion m ~ 1.5 results in t^2 — 700 yr. Concerning zone 1, 
we note that along with the compact radio structures with 
declining brightness observations show also a large num- 
ber of brightening components, and the timcscalcs of the 
flux variations in the individual knots can be as small as 
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few tens of years (Tuffs 1986| , AR96). An estimate of the 
average time fsi of the magnetic field decrease in CBSR 
components is thus problematic. Therefore we consider £bi 
as a free model parameter, and neglect adiabatic losses of 
electrons there (if any - especially as they are supposed 
to be the sites of effective acceleration). 




Fig. 7. The rate of secular decline of the total flux of 
synchrotron radiation shown in Fig. 4, calculated for char- 
acteristic timescales (see text) £32 = 650 yr and £ a d2 = 
1000 yr for zone 2, and different timescales for the mag- 
netic field decline in zone 1: t&i = 115 yr (dot-dashed), 
130 yr (solid), 150 yr (dashed). The diam ond shows the 
recent result of O'Sullivan & Gr een (|l999|) at 15 GHz, all 
other data points are from Rees ( 1990 ). The heavy dashed 
straight line corresponds to the earlier 'best fit' estimate 
(see text). 

Calculations of the secular decline rate of the fluxes at 
different frequencies are presented in Fig. 7. The straight 
heavy dashed line shown corresponds to the decline rate 
deduced by Dent et al. ( |1974[ ) and Baars et al. ( 1977 ) from 
comparison of the decline rate ~ 1.29 %/yr measured then 
at 81 MHz by Scott et al. fll969j ) with the data in the 1- 
10 GHz region. Later measurements have shown that the 
decline rate around v ~ 100 MHz is at the same level 
~ (0.8 — 1) %/yr as at v ~ 1 GHz, which has placed doubt 
on the frequency dependence of the decline rate in general 
(Rees [l990 Hook et al. 1992). However, the recent result 
of O'Sullivan & Green ( 1999] ) that the flux decline rate at 
15 GHz is about 0.6%/ yr, seems to confirm the effect of 
slower decline of radio fluxes at high frequencies. 

Figure 7 shows that practically all the observational 
data can be well explained assuming the average magnetic 
field decline time in zone 1 of order t-Qi ~ 120 — 150 yr. 
Note that a significant drop of the secular decline rate 
predicted at v < 30 MHz is connected with the gradual 
decrease of synchrotron self-absorption in zone 1. 

5. CBSR components as acceleration-passive sites 

Consider now the hypothesis that acceleration of electrons 
takes place only in zone 2 and is negligible in zone 1. In this 
case the source function Q2(E) for zone 2 is in the form 



of Eq.(18), while Qi(E) for zone 1 is defined only by the 
rate Q2 1 (E) of advection of relativistic electrons from the 
plateau region. The required for zone 1 relativistic electron 
density enhancement could be then supposed as being due 
to adiabatic compression of the zone 2 electrons on the 
front surface SL of fastly moving CBSR components, and 
their subsequent accumulation in zone 1 over timescales 
At ~ r osc . Neglecting in Eq.(7) the diffusion term on the 
surface SL, in compliance with the assumption of neg- 
ligible diffusive acceleration there, and retaining only the 
convection and the adiabatic compression terms, the equa- 
tion describing the transport of the electrons through 5_ 
reads: 



dfp _ dux_ pdfp 
dx dx 3 dp 



= 



(22) 



where f p = f p (x,p,t) is the momentum distribution func- 
tion, x is the spatial coordinate and u x is the plasma 
speed along the x-axis perpendicular to the compression 
layer. General solution to this equation reads: f P (p,x) = 
Foipiux/uo) 1 / 3 ], where F (p) = f p (p,x ) is an arbitrary 
initial function. Connecting the electron distribution func- 
tions across the layer, the injection rate of relativistic elec- 
trons from zone 2 into zone 1 is found: 



Q 21 (E) ~ N 2 (E p' 1 ' 3 ) p- 1 ' 3 S- u 2 Vi x , 



(23) 



where p c = U2/U1 is the mean compression ratio, and u 2 
and Mi correspond to the plasma speeds of the plasma in 
the upstream and downstream regions in the rest frame of 
the CBSR components. 

Considering possible values of 112, note that for the 
mean radial expansion age of radio knots t cxp = r/v r 
withi n (550 - 1000) yr (Tuffs |1986j ; Anderson & Rudnick 
1995) the characteristic speed of the knots may be es- 
timated as «kn — (1700 — 3200) km/s. Taking into ac- 
count that surrounding plasma in the shell is moving in 
the same direction, and as evidenced by a shorter X-ray 
expansion t ime scales ~ 600 yr (Koralesky et al. 1998 ; 
Vink et al. 1998] ), appears to be overtaking most of the 
compact radio knots, the value 112 ~ 2000 km/s seems a 
reasonable estimate for the relative speed of radio knots 
with respect to the surrounding medium. For the bright 
radio ring i cxp ~ 950 yr (Tuffs [1986 ), corresponding to 
^ring — 1700 km/s, while the speed of freely expanding 
ejecta upstream of the reverse shock, i.e. at r < i? r i ng , 
is about (5000 — 5300) km/s. Thus, for the radio ring 
U2 ~ 3500 km/s, and we can take 3000 km/s as a rea- 
sonable maximum value of 112 for the zone 1. 

The compression ratio p c can be approximated as the 
ratio of magnetic fields p c ~ B\jB2- The energy distribu- 
tions of the radio electrons can be connected as 



Ni{E*) Yl 
N 2 (E*) ~ V 2 



£2-1 

Bi\ 3 U 2 T* 



2i?i 



(24) 
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using the relation Ni(E*) ~Q 2 i x r*. Eqs. (3) and (24) 
result in 



Bi 
B 2 



5_g2 + l 



Ji 2R 1 V 2 (1Bi\ 2 



(25) 



where A/3 = /3j — (3 2 . Since the volume filling factor of 
the compact zone 1 structures is Vi/V 2 ~ 0.5% or less, 
and since for effective modification of the electron spectra 
by escape one needs < 0.1 to — 30 yr, the ratio of the 
magnetic fields in the two zones should be rather high, 
Bx/Bi > 10. 

The results of numerical calculations are shown in 
Fig. 8. The total flux seems to fit the data practically 
equally well as in Fig. 4 earlier. However, a closer look 
to these figures reveals significant differences. First of all, 
in Fig. 8 the agreement with the radio data is reached due 
to a much higher contribution of zone 1 to the overall 
flux. In particular, at 5 GHz the partition of the fluxes 
corresponds to J\/ J 2 ~ 3, instead of J\/J 2 — 1- Secondly, 
although here we have assumed high 5 = 0.75, the calcu- 
lations show that the spectral index of the electrons swept 
up from zone 2, j3 2 — 2.2, was steepened in zone 1 only 
to values j3\ < 2.6, well below the 'expected' value 2.95. 
Obviously, this is far insufficient for explanation of the 
observed steep spectral indices of the radio knots. 




Fig. 8. Radiation fluxes expected in the scenario when ac- 
celeration occurs only in the zone 2 (diffuse plateau), but 
not in the compact zone 1. The solid lines show the total 
fluxes, while the dashed and dot-dashed lines correspond 
to the contributions from the zone 1 and zone 2, respec- 
tively. Model parameters are: Bi = 3.8 mG, B 2 = 0.2 mG 
/3 acc = 2.2, E c = 30,TeV, r ast = 20 yr, 8 = 0.7, K = 200, 
Ri = 0.06 pc, and u 2 = 3000 km/s. 

The reason for difficulties arising in the scenario that 
assumes only passive compression, but not active accel- 
eration of electrons in zone 1, can be understood from 
Fig. 9 where we show the spatial densities rii t2 (E) of the 
electrons in zones 1 and 2. The dot-dashed line shows the 
spectrum of the electrons which would be formed in the 
zone 1 if one neglects the term oc n 2 /Tdif in Eq.(14). This 
spectrum is steep. But when the effect of spatial density 




Fig. 9. The energy distributions of the electrons formed in 
zone 1 (solid line) and zone 2 (dashed line) in the frame- 
work of 'adiabatic compression' scenario for zone 1 elec- 
trons, calculated for parameters in Fig. 8. The dot-dashed 
curve corresponds to the the energy distribution of elec- 
trons in zone 1 when the return diffusive flux of electrons 
from zone 2 into zone 1 is neglected. 

gradients between the two zones is taken into account cor- 
rectly, the situation dramatically changes. The energy de- 
pendent propagation can result in a significant steepening 
of the source spectra only in the regions of higher spa- 
tial densities of relativistic particles. Otherwise diffusive 
exchange of particles tends to equalize ni(E) with n 2 (E) 
resulting in /3i — ► (i 2 independently of S, until the radiative 
losses become important. 

Thus, for an effective modifications of the initial source 
spectrum of the radio electrons the condition n\{E) 3> 
n 2 (E) should be satisfied. The ratio of electron densities 
expected at 1 GeV can be estimated from Eqs. (25) and 
(26) (assuming (3 2 ~ 2.2) as: 



«*2 



0.2 



0.2 



U 2 T» 



/2*\ 
VlmG J 



-0.1 A/3 



.(26) 



For the parameters used in Fig. 8 this equation predicts 
n*i/n*2 ~ 2, in agreement with numerical calculations, 
which is smaller by a factor of 15 than similar ratio in 
Fig. 3. Calculations show that it is not easy to increase this 
ratio significantly, remaining within the 'adiabatic com- 
pression' scenario. 

Thus, a self-consistent explanation of the observed ra- 
dio fluxes of Cas A is very problematic, unless we assume 
an effective acceleration of electrons in the CBSR compo- 
nents in order to build up sizeable gradients in the spatial 
distribution of the radio electrons in the source. 

6. Results 

The spatially non-uniform source model, even in its sim- 
plest two-zone form, allows us to unify into a single picture 
many observational data on the broad-band nonthermal 
radiation of Cas A. The following is the summary of ob- 
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servational features of Cas A which could be explained in 
the framework of this model. 

The brighter radio structures tend to be steeper (AR96) 
to the extent that the energy density of relativistic elec- 
trons there would be higher. In principle, even assuming 
the same hard power law index for the acceleration, e.g 
/3 acc ~ 2.2, and S ~ 0.6 for the escape of the electrons, it 
is possible to explain the observed variations of the spec- 
tral indices of individual compact radio structures from 
a 2 ~ 0.6 to ai ~ 0.9. Such variations can be connected 
with two effects. 

(a) The intrinsic index of some knots will be flatter 
than the maximum possible ct\ — a 2 + 5/2 if the local 
gradient n\jn 2 is not sufficiently high, and if the charac- 
teristic escape time r» of GeV electrons from these knots 
is larger than ~ 30 yr. The impact of the latter effect can 
be seen in Figs. 10 and 11, where we show the fluxes and 
spectral indices calculated in the modified two-zone model 
approach, subdividing the CBSR structures in two groups 
and assuming 2 different spatial/temporal scales for larger 
('ring') and smaller ('knot') components in zone 1. 
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Fig. 11. The spectral indices of the radiation components 
shown in Fig. 10 

(b) The intrinsic spectrum of individual knots could be 
steep, up to ai ~ 0.9, but the observed spectrum may be 
flatter, depending on the degree of flux contamination by 
the flat-spectrum radiation from the diffuse plateau along 
the line of sight to the knot. This effect may be relevant 
mainly to those cases when the background subtraction 
is problematic (e.g., for the structures with insufficiently 
high contrast to the plateau). 

The strong correlation between the spectral index and the 
projected position of radio knots observed by AR96 could 
be connected with both effects (a) and (b). Indeed, the 
brightness of the diffuse plateau seems to decrease from 
the radio ring towards the outer edge of the shell more 
rapidly than one would expect merely on the base of the 
projection effect in a spherical geometry (see Fig. 7 from 
ARLPB). This implies that the density n 2 of radio elec- 



trons could be higher at distances closer to the ring. Then 
the steepest knots would be found predominantly closer to 
the edge of the plateau, because (i) the local ratio n\ /n 2 
might be higher far from the reverse shock, and (ii) contri- 
bution of the flat spectrum plateau emission to the steep 
spectrum flux of the knots is smaller in directions to the 
projected periphery of the shell. 

The cutoff of the radio spectrum below 20 MHz is ex- 
plained by the synchrotron self-absorption of the flux of 
zone 1. This interpretation is possible basically because 
the model predicts that the intrinsic fluxes of the com- 
pact bright structures are much steeper than the spectra 
in the diffuse plateau. Then at low frequencies the flux of 
CBSR components should dominate the total radio emis- 
sion. Therefore synchrotron self-absorption of only this 
flux, which becomes possible due to high density of parti- 
cles (and fields) in those compact structures, is sufficient 
to comply with the data. 

It is worth noticing that due to variations of the phys- 
ical parameters in the individual CBSR components con- 
tributing to zone 1 flux, one could generally expect that 
the synchrotron absorption of that flux would be in re- 
ality smoother than it is shown in Fig. 4 where we have 
assumed the same parameters for all counterparts of zone 
1. The position of the characteristic cutoff frequency for 
an individual component j is found from the condition 
Tj(y) — 1 where Tj — RjK(v) is its opacity. From Eqs. 
(3) and (4) it follows that Tj(u) cx A j iy-^ 2+l3 ^ 2 \ with 
Aj = RJ 2 fjBj 5 where fj is the knot flux (luminosity) 
at some fixed frequency (say 1 GHz). Because of the very 
strong dependence of tj on v, a significant broadening of 
the synchrotron turnover frequency for the ensemble of 
CBSR components (say by a factor 2) would be expected 
only if the luminosity-weighted dispersion a of the param- 
eter Aj is very large (a ~ 10). The fact that the observed 
turnover is sharp may then imply a quite reasonable possi- 
bility that the distribution of the parameters {Aj } is not 
far from 'Gaussian noise'. In that case one would normally 
expect a ~ 1 (or perhaps even less, given the expected cor- 
relation of fj and Rj), so the synchrotron turnover posi- 
tion would be effectively dispersed, or broadened, only by 
a factor roughly (1 + a) - 3 ~ 1.2, i.e. about ±20% around 
the mean position. Note that the overall flux in Fig. 10, 
where 2 significantly different sizes for the large and small 
components are assumed, show practically the same sharp 
cutoff below 20 MHz as in Fig. 4. 

Observations at low frequencies by KPDE, at a resolution 
of about 20", show that the spectral index of Cas A be- 
tween 333 MHz and 1.38 GHz is approximately constant, 
a ~ 0.75, at angular radii 9 < 100", but that it is quickly 
increasing to a < 0.95 as 8 — > 150". In principle, this 
effect might be connected with the steep spectrum radio 
knots which have been detected by AR96 but cannot be 
resolved in the low frequency maps of KPDE. Because of 
the rapid decline of the plateau brightness at large angu- 
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Fig. 10. Fluxes expected in the framework of the modified two zone model, when in zone 1 more compact radio knots 
(dashed line) and larger fragments of the radio ring (thin solid line) are separated. The 3-dot-dashed line shows the 
total flux from zone 1, and the dot-dashed line shows the flux from zone 2. The total synchrotron flux is shown by 
the heavy solid line. The model parameters are the following. For the knots: i?i.kn = 0.033 pc, r* = 15 yr, 5 = 0.6, 
iikn = 180; for the ring fragments i?i. r i ng = 0.1 pc, r* = 35 yr, S — 0.6, i^ r i ng = 20. The same /3 acc = 2.2 and cutoff 
energy E c — 35 TeV for the accelerated electrons, and magnetic fields £?i.kn = B\ y mg = 1.5 mG in the knots and the 
ring is assumed. The magnetic field B2 = 0.36 mG. 



lar distances from the ring, the CBSR components (the 
radio knots and bow shocks) increasingly contribute to 
the 'diffuse' overal flux closer to the periphery. We do not 
exclude that perhaps in the framework of a more spatially 
structured model (that would allow large scale inhomo- 
geneities in the plateau region itself) one could expect 
also some steepening of the intrinsic diffuse flux. However, 
given much slower and less efficient spectral modifications 
for the large scale inhomogeneities, the propagation ef- 
fects alone would be able to explain only rather moderate 
steepening of the intrinsic plateau emission. If the accel- 
eration is efficient so that the primary spectrum is really 
hard, the steepening to values a ~ 0.9 seems to require a 
dominating contribution of the compact structures in the 
overall flux detected by KPDE. 

The two-zone model suggests a possible interpreta- 
tion for the 'discrepancy' between the spectral index a = 
(0.7-0.75) found by KPDE at the radius of the bright ra- 
dio ring, and observations of WDS90 who found a ~ 0.83 
if only the bright components at approximately the same 
radial distances are considered. Both results might be un- 
derstood if we take into account that the spectral index 
found by KPDE corresponds to the flux of two zones inte- 
grated along the geometrical circles of different fixed radii, 
and hence would be more affected by the flat plateau flux 
than the result found in WDS90 (compare the heavy and 
thin solid lines in Fig. 11). 

The flattening of the radiation spectrum observed at mil- 
limeter wavelengths (Mezger et al. 1986) is a natural con- 



sequence of the 'flat + steep' representation of the total 
flux in the two-component approach, with approximately 
equal contributions from these components around 5 GHz. 
The two-zone model also predicts that the flux measured 
by Tuffs et al. (1997) around 6 \xm should be predomi- 
nantly synchrotron in origin. 

The unusual dependence of the secular decline of ra- 
dio fluxes on frequency, which is at a constant level at 
v ~ 40 MHz — 1 GHz but seems dropping at higher fre- 
quencies, is explained by varying contributions of the flat 
and steep spectrum components at different frequencies. 
This interpretation suggests that on average the net flux 
of the bright compact components in Cas A drops faster 
than the diffuse plateau emission, in agreement with ob- 



servations of Tuffs (1983) 



The synchrotron origin of the X-rays above 10 keV can be 
explained if electrons in the compact CBSR components 
are accelerated to energies of few tens of TeV. These en- 
ergies can be reached, e.g., by diffusive shock acceleration 
in the Bohm limit, which seems seems a plausible mech- 
anism for production of relativistic electrons at the re- 
verse shock presumably connected with the bright radio 
'ring'. In principle, for the radio knots as well the diffu- 
sive acceleration at the bow shocks could result in the high 
electron energies needed for X-ray emission. Perhaps, an- 
other possibility for acceleration of electrons in the radio 
knots could be connected with reconnection of the mag- 
netic field lines strongly amplified, as shown by Jones et 
al. (1994), in thin turbulent layers behind the bow shocks 



16 



A. M. Atoyan et al.: Energy-dependent propagation effects in Cas A 



of dense gas 'bullets' at the stage of their fast disruption 
by Rayleigh- Taylor and Kelvin-Helmholtz instabilities. 

Hard spectral indices for accelerated particles weaken the 
constraints on the energetics of relativistic electrons in 
Cas A. For example, in the case of /3 acc = 2.24, and B\ = 
1.5 x 1CT 3 G and B 2 = 4x 1CT 4 G used in Fig.4, the total 
energy W c = 2.7 x 10 48 erg. The energy W e needed in 
the case of (3q — 2.54 for a uniform shell with the same 
magnetic field Bq = B2 is by a factor of 10 larger. 

7. Conclusions 

The radio fluxes of the prototype young SNR Cas A show 
a large variety of spectral and temporal features which 
are very difficult to incorporate into one self-consistent 
picture, if one remains in the framework of a single-zone 
(i.e. uniform) model approach and thus has to attribute 
the spectral indices deduced from radio observations to 
the source spectra of accelerated electrons, since radiation 
losses cannot modify the energy distribution of radio elec- 
trons during the lifetime of Cas A. The interpretation of 
the radio data essentially changes if we take into account 
that energy dependent propagation of relativistic particles 
is able to modify their energy distribution in a spatially 
inhomogencous source on timescales much shorter than 
the radiative loss time. The efficiency and extent of these 
modifications depend on the strength of the gradients in 
the spatial distribution of the particles in the source. 

The simplest spatially inhomogeneous model is the 
one where the radio source is subdivided into two zones 
with different energy distributions of relativistic electrons 
and different magnetic fields. The basic assumption of the 
model is that the number density of radio electrons in 
the compact bright radio components (predominantly in- 
cluded in zone 1) is much higher than in the surround- 
ing diffuse shell between the reverse shock and the blast 
wave (zone 2). The fulfillment of this assumption is contin- 
gent on an efficient acceleration of electrons in these radio 
bright components. In principle, acceleration of electrons 
in the diffuse shell is not needed at all. Note, however, that 
the model does not exclude a contribution from such accel- 
eration provided that it is not so powerful as to noticeably 
reduce the gradients in the spatial density of electrons es- 
tablished by acceleration in the CBSR components, which 
would otherwise diminish the efficiency of spectral modi- 
fications in zone 1. 

The possibility for exchange of particles between these 
zones on an energy dependent timescale t csc oc E~ 5 , 
with 5 ~ 0.5 — 0.6, allows us to suggest a single hard 
power law injection spectrum of electrons with an index 
/3 acc ~ 2.2 — 2.3, implying an efficient acceleration pro- 
cess. Then, depending on the relativistic electron density 
contrast and physical parameters of individual component 
j, the leakage of electrons from those CBSR components 
can result in the steepening of the injection spectrum up 
to (3ij < 2.9, in agreement with the maximal spectral in- 



dices of the knots ctij < 0.95. The energy distribution of 
radio electrons N 2 (E) in the extended plateau will show 
the hard power-law index of injection, /? 2 ~ /3 aC c- As a gen- 
eral remark we note that energy-dependent propagation in 
any non-uniform medium would always tend to flatten the 
energy distribution of particles in the regions with lower 
density, i.e. which are typically more extended. 

The model leads to a number of predictions. 

The two-component decomposition of the overall flux, 
with Ji ~ J 2 at 5 GHz, predicts that the sites of current 
or very recent acceleration will become more pronounced 
in the high resolution (few arcsec) maps at lower frequen- 
cies, where the contribution of the diffuse plateau to the 
overall flux will decrease. This should result in a signifi- 
cant increase of the brightness contrast for the remaining 
compact structures (which belong to zone 1), reaching its 
maximum around 40 MHz (i.e. at the minimum frequen- 
cies not affected by synchrotron self-absorption). 

At frequencies < 30 MHz the bright compact struc- 
tures will quickly disappear, but the plateau emission, in 
the form of a weak diffuse shell, may become dominant 
again. This shell would probably have an apparent ra- 
dius less than 150" if the emissivity of the shell is indeed 
strongly decreasing towards the blast wave, implying very 
low brightness at large radii. 

Correlations between the spectral index and the geo- 
metrical thickness of the radio knots of similar brightness 
could be expected on the high resolution maps of Cas A 
taken at 1.4-5 GHz. The radio knots with a smaller thick- 
ness, but the same brightness, would tend to be steeper 
since the escape times in those knots could be smaller. 
Note however that the escape time is not the only param- 
eter affecting the efficiency of spectral modifications. In 
particular, the 'steepness-compactness' correlation could 
be different for knots in the fading stage (which are 'older', 
with declining injection of new particles) and in the bright- 
ening stage. Besides, such a correlation could be signifi- 
cantly affected by different magnetic fields in the knots. 
Nevertheless, the search for correlations between spectral 
index and compactness of the knots seems worthwhile. 

At frequencies above 10 GHz the brightness contrast 
will decrease, and the CBSR structures may appear less 
pronounced. The total emission will be dominated by the 
diffuse flat-spectrum plateau, and the spectral index may 
decrease to a « 0.65 (see Fig. 11). The structure of zone 
2, in particular the outer edge of the shell, will be better 
discernible. 

An important prediction concerns the character of the 
long term evolution of the radio pattern of Cas A. One can 
expect that during the next At ~ £bi ~ (150 — 200) yr, 
when the total emission will become dominated by the 
diffuse plateau and the compact radio structures may be- 
come less common, the remnant will have a much more 
uniform radio shell with a spectral index ~ 0.6, and Cas A 
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may become similar to Tycho's SNRQ (e.g. see Klein et al. 



1979). Interestingly, the overall efficiency of electron accel- 



eration at that time might be even lower than at present, 
while the radio spectrum will be significantly harder re- 
flecting the spectrum of the acceleration which is tak- 
ing place presently. More generally, the energy dependent 
propagation of radio electrons in a spatially inhomoge- 
neous medium can explain the trend (see e.g. Green 1988; 
Jones et al. 1998) that young clumpy shell-type SNRs 
often exhibit radio spectra that are significantly steeper 
than those of older ones showing a typical spectral index 
a ~ 0.5. 

At X-ray frequencies the appearance of Cas A may 
be very different in the soft X-ray and hard X-ray do- 
mains. At photon energies below 10 keV the radiation is 
dominated by thermal emission of the gas in the extended 
region from the reverse shock up to the blast wave. Apart 
from the thermal diffuse emission we could expect a no- 
ticeable contribution of the nonthermal X-rays from com- 
pact radio structures. This radiation could be more eas- 
ily distinguished from the thermal emission in the case of 
compact bright radio knots at the periphery of Cas A, if 
the X-ray spectra would appear relatively flat in the 1- 
5keV region (see Figs. 4 and 11) and the line emission 
would be deficient. These observations will be possible 
with the high spectral and angular resolutions of XMM 
and ASTRO-E. 

In hard X-rays above lOkeV the appearance of Cas A 
may significantly change. In the case of a nonthermal ori- 
gin of this radiation, all acceleration sites of the highest 
energy electrons may become clearly visible. Note, how- 
ever, that the angular resolution of the detectors needed to 
reveal significant spatial changes at these energies must be 
better than 10". The observation of such hard nonthermal 
X-ray fluxes from the compact radio knots and radio ring 
would impose an upper limit on the magnetic field in these 
structures B\ < 3 niG, with the probable value expected 
in the range B\ ~ (1 — 2) mG. The model predictions for 
magnetic fields in the diffuse shell are B2 — (0.3 — 0.5) mG. 

At energies above 100 keV we predict very flat spectra 
of the soft 7-ray fluxes due to bremsstrahlung which per- 
haps would be observable for the forthcoming ASTRO-E 
and INTEGRAL telescopes. 

For the high energy 7-rays, E > 100 MeV, fluxes at 
a level of 0.1-0.5 of the EGRET flux upper limits are 
expected. These fluxes should be easily detected by the 
GLAST instrument. 



3 An alternative explanation for the contrasting brightness, 
spectral and morphological characteristics of Tycho's SNR 
could be that, as the remnant of a type la event, this source 
never contained compact efficient accelerators as appear to be 
now present in the ejecta - circumstellar interaction of Cas A, 
and that all particle acceleration in Tycho's SNR has occurred 
at the blast wave. 



In summary, the steepness of radio spectra of bright and 
compact structures in clumpy radio sources is not an indi- 
cation of inefficient acceleration, but rather a natural con- 
sequence of very efficient acceleration which builds up high 
spatial gradients of relativistic electrons in those sources 
increasing the efficiency of spectral modifications due to 
their energy dependent escape. 
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